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ABSTRACT 

We report measurements of anisotropy in the cosmic microwave background radiation over the multi- 
pole range i ~ 200 — > 3500 with the Cosmic Background Imager based on deep observations of three 
fields. These results confirm the drop in power with increasing t first reported in earlier measurements 
with this instrument, and extend the observations of this decline in power out to i ~ 2000. The de- 
cline in power is consistent with the predicted damping of primary anisotropics. At larger multipolcs, 
£ = 2000-3500, the power is 3.1c greater than standard models for intrinsic microwave background 
anisotropy in this multipole range, and 3.5(7 greater than zero. This excess power is not consistent 
with expected levels of residual radio source contamination but, for er 8 > 1, is consistent with pre- 
dicted levels due to a secondary Sunyaev-Zeldovich anisotropy. Further observations are necessary to 
confirm the level of this excess and, if confirmed, determine its origin. 
Subject headings: cosmic microwave background — cosmology: observations 



1. introduction 

Sub-horizon scale fluctuations in the Cosmic Mi- 
crowave Background (CMB) provide a direct view of 
simple, causal physical processes in the early universe. 
In standard cosmological models, the dominant pro- 
cesses are acoustic oscillations of the primordial plasma 
and photon-diffusive damping, which give rise to a har- 
monic series of peaks in the CMB anisotropy spectrum 
mo dul ated by an exponential cutoff on small angular 
scales llBond fc Efstathioull987tlPeebles fe Yulll970tlsilkl 
mvaev fc Zel'dovichlll970ri . Measurements of the 
CMB power spectrum in this regime provide strong con- 

1 Current Address: National Radio Astronomy Observatory, 
P.O. Box 2, Green Bank, WV 24944 



straints on cosmological parameters ijWhite et alJ ll9941. 
determine the nature and in itial conditions of the fluc- 
tuations IHu fc White! llQQfij), and provide fundamen- 
tal te sts of particle physics ijKamionkowski fc Kosowskil 

1999) . A number of experim ents have detected peaks 
in the anisotropy spe c trum l|4 ^ Bernarcjis et al.| I2QQQJ: 
I Halverson et al.l l2002t lHananv et all 120001: ILee et alJ 
boOll iMiller et all 119991: iNetterfield et all 12002(1 . which 
constrain the cosmological models through measure- 
ments of the first, second, and possibly the third acous- 
tic peaks. Observations at high multipoles (t ~ 500 — > 

2000) , where the physics is strongly affected by photon 
diffusion and the thickness of the last scattering region, 
provide independent constraints on these fundamental 
parameters. At even higher multipoles (£ > 2000), 
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secon dary effects such as the Sun yaev-Zeldovich Effect 
(SZE; Sjrnv^iev are expected to dom- 

inate <|Cole fc Kaiserlll988HReDhaelilll981|) and hence of- 
fer the prospect of studying the formation of large-scale 
structure at recent times. 

This paper is one in a series reporting results from the 
Cosmic Background Imager (CBI). We have previously 
reported results between £ — 600 and £ = 1 200 in which 
we de tected a damping tail in the spectrum ijPadin et al.l 
120011 hereafter Paper I). These measurements support 
the standard theoretical models of CMB anisotropics. In 
the present paper, Paper II, we present measurements 
with the CBI taken from January through December of 
2000. These observations extend our determination of 
the CMB anisotropy spectrum out to £ = 3500. The re- 
sults are derived from deep integrations on three pairs of 
pointings. We present a brief overview of the instrument 
and site in § [21 and discuss our observing technique in 
§ The data analysis methods are presented in § 0J in- 
cluding, in S 14.21 a discussion of the maximum- likelihood 
analysis used to determine the power spectrum and, in 
§ 14.31 a discussion of the method of removing the dis- 
crete source foreground. We present our results in § 
and review our conclusions in § . 

The remaining pap ers in the series cove r the CBI mo- 
saic power spectrum ({Pearson et al.ll2003L hereafter Pa- 
per III); t he implementation of the maximum likelihood 
analysis l)Mvers et al.l [2003, hereafter Pa per IV); the 
cosmological interpretation of our results l)Sievers et alJ 
I2003L hereafter Paper V); and a possi ble interpretation 
of the excess power observed at high-£ ijBond et al.ll2003l 
hereafter Paper VI). 

2. THE INSTRUMENT AND SITE 

A detailed description of the CBI can be found in 
the paper bv iPadin et al.l 1(2002), so we summarize here 
only important aspects of the instrument design. The 
CBI is a planar synthesis array of 13 0.9-meter diame- 
ter Cassegrain antennas mounted on a 6-meter diame- 
ter tracking platform. Short baselines arc susceptible to 
contamination from cross-talk, so the antennas are sur- 
rounded by cylindrical shields which provide more than 
110 dB of isolation. Further rejection of cross-talk is ac- 
complished through the use of differenced observations 
(see © ■ A key feature of the CBI is its sensitive broad- 
band receivers, which are based on indium phosphide 
high electron mobility transistor (HEMT) amplifiers op- 
erating in the 26-36 GHz band. The system tempera- 
tures measured on the telescope, including CMB, ground, 
and atmosphere, are typically ~ 30 K. The receiver out- 
puts are combined in an analog correlator with 10 1-GHz 
bands. A 1 kHz phase switching scheme is used to reject 
cross-talk in the electronics. 

The antenna platform is on a 3-axis mount with az- 
imuth, elevation, and parallactic angle axes; the eleva- 
tion is restricted to > 42°. The parallactic angle rota- 
tion provides the ability to track the rotation of the sky 
throughout the observation so that the (u, v) orientation 
of a baseline on the sky is fixed, where u and v are the or- 
thogonal components of the baseline length measured in 
wavelengths. The platform is moved through additional 
discrete steps in orientation to increase the (u, v) cover- 
age, test for false signals, and increase ground rejection. 
The antennas can be placed in different locations on the 



platform, allowing the configuration to be matched to a 
range of science goals. 

The CBI is located at an altitude of 5080 m in 
the Atacama desert of northern Chile, near the pro- 
posed site for the Atacam a Large Millimeter Array 
( Radf ord & Holdawavl [19981. This high dry site was 
chosen in order to reduce atmospheric brightness fluc- 
tuations which would otherwise limit the sensitivity of 
the CBI. During the calendar year 2000, ~ 35% of the 
nights were lost to weather and ~ 5% to equipment 
malfunction. Observing conditions were best from mid- 
August through mid-December, when no nights were lost 
to weather. In the nights when observing conditions were 
good, there is no evidence of atmospheric effects in the 
CBI data except for rare instances which typically cor- 
relate with the appearance of visible, low-lying clouds. 

3. OBSERVATIONS 

The data for the results presented here were obtained 
from 2000 January 11 to December 12. During this pe- 
riod both deep observations of individual fields and mo 
saiced observations of multiple pointings were made. We 
present in this paper the results for three deep fields, 
two of which are part of the mosaics; the CBI mosaic 
data, images, and power spectra are discussed in Paper 
III. We refer to these three deep fields collectively as the 
CBI deep fields, and individually as deep 08h, deep 14h, 
and deep 20h. We obtained 42, 9, and 49 good nights of 
observing on these fields. Thus most of our results come 
from the deep 08h and deep 20h fields, with the deep 14h 
providing a somewhat weaker independent check. All of 
the deep fields were chosen to have low point-source con- 
taminat ion levels using the N RAO Very Large Array Sky 
Survey ijCondon et al.ll1998l NVSS). 

In order to eliminate the influence of solar radiation in 
distant sidelobes, observations were made only at night. 
For similar reasons, no observations were used for which 
the field was less than 60° from the moon. Our results 
are insensitive to the precise lunar cut beyond 60°. As 
discussed in Paper I, the CBI has no ground shield and 
observations are differenced in order to remove signals 
due to ground pickup on the short baselines. To accom- 
plish this, two fields, designated lead and trail, separated 
by 8 minutes in right ascension, are tracked across al- 
most identical ranges of azimuth and elevation. Ground 
pickup and other contaminating signals are cancelled in 
the lead — trail difference as long as they are stable over 
an 8 minute time-span. The potential residual of this 
in the differenced data has been shown in Paper I to be 
< 1.3% of the cosmic signal. 

The coordinates of the lead and trail fields for the 3 
CBI deep fields, the total integration times on lead and 
trail combined, and the Galactic latitudes are shown in 
Tabled The integration times shown are computed after 
all data filters have been applied. 

For the observations reported in this paper, four an- 
tenna configurations were used. The first configuration, 
on which the results of Paper I were based, was a ring 
configuration, with most of the receivers near the perime- 
ter of the platform. This provided fairly uniform (u, v) 
coverage, many long baselines for robust calibration, and 
good access to the antennas and instrumentation on the 
platform during the testing which accompanied initial 
observations. The second and third configurations were 
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Name R,.A.j ea ^ (J2000) R.A. trail (J2000) Dccl. (J2000) b t int 



deep 08h 08:44:40 08:52:40 -03:10:00 24° 131.2 h 

deep 14h 14:42:00 14:50:00 -03:50:00 40° 24.0 h 

deep 20h 20:48:40 20:56:40 -03:30:00 -28° 115. 7h 



I I 08 h 14 h I I 20 h 




Multipole I 



Fig. 1. — The distributions of the weights of the visibilities as 
a function of i for the three deep fields. Since the CBI sensitivity 
does not vary much with baseline or frequency, these weights are 
roughly proportional to the on sky integration time as a function 
of L 



better optimized for the CBI mosaic observations. These 
more compact configurations provided higher sensitivity 
in the I = 600-1500 range, and were well-suited to low- 
rcdshift Sunyaev-Zeldovich observations. In addition, 
the third configuration provided nearly redundant base- 
lines at different frequencies which enabled us to make 
a more accurate determination of the radio frequency 
spectrum of the signals reported in Paper I. The last 
configuration was an extended configuration with more 
uniform coverage than the first, more long (> 5 m) base- 
lines, which are useful for point source monitoring, and 
greater sensitivity for the frequency spectrum determina- 
tion out to £ = 1400. The distributions of the weights of 
the data as a function of £ is shown in Figure The 08h 
field was observed in the first, second and fourth config- 
urations; the 14h field in the first configuration; and the 
20h field in the third configuration. 

In the first three configurations, all but one of the re- 
ceivers were configured to be sensitive to left circular 
polarization (LCP); the remaining receiver was sensitive 
to right circular polarization (RCP) . In the last configu- 
ration all antennas were configured for LCP. 

4. DATA ANALYSIS 

4.1. Calibration and Editing 

The equatorial sky coverage and relatively large an- 
tennas of the CBI enable us to use a variety of ce- 
lestial objects as calibrators, including planets, super- 
nova remnants, radio galaxies and quasars. This is con- 
venient for comparing the CBI flux density scale with 
other measurements (see below) and for accurately cal- 
ibrating the data. CBI observations are calibrated us- 
ing nightly measurements of primary calibrators, cho- 



sen for their strength and lack of variability at these 
frequencies (Jupiter, Saturn, Tau A [3C144, the Crab 
Nebula], and Vir A [3C274]), and a set of secondary cal- 
ibrators (3C279, 3C273, J1743-038, B1830-210, and 
J1924— 292). Our estimates of the primary calibra- 
tor flux densities are based on CBI measurements of 
these sources relative to Jupiter assuming an apparent 
Ravleigh- Jea ns brightness temp erature Tj up = 152 ± 5 K 
at 32 GHz i|Mason et al J 119991. This value is uncor- 
rected for the occultation of the CMB by the planet, and 
applying this correction yields a Rayleigh- Jeans bright- 
ness temperature of 154 ± 5 K. Thi s is also the basis o f 
the Very Small Array calibration (Wat son et al J 12002'!. 
Since Jupiter is known to have a non-thermal spectrum 
at these frequencies, and the precise spectral shape is not 
independently well-determined, this flux density scale is 
transferred to each of the CBI bands by observations of 
Tau A, for which the iBaars et alJ l)1977l) spectral index 
a = —0.299 is assumed, where S oc v a . There is good ev- 
idence th at the power law ind ex extends to high frequen- 
cies (e.g.. iMezger et aj]ll986D . The flux densities of the 
secondary calibrators are bootstrapped from the nearest 
primary calibrator observations, and are used when the 
primary calibrators arc cither not visible or are too close 
to the moon. Observations of the primary calibrators 
over the year 2000 indicate random calibration errors of 
~ 3% night-to-night. The residuals of a polynomial fit 
to the light curves of the secondary calibrators show a 
similar scatter. I n light of this and th e 3.3% calibration 
uncertainty of the lMason et al.l Jl999) value for the tem- 
perature of Jupiter we assign a calibration uncertainty of 
5% to our data. This corresponds to a 10% uncertainty in 
bandpower (/iK 2 ). These include the beam uncertainties 
discussed below. 

We have two comparisons of our flux density scale with 
that of the NRAO Very Large Array 2 (VLA) at 22.4 
GHz. First, we compare CBI observations of Mars with 
the m odel for the Mars temperature which is used at the 
VLA (|R.udvlll987tl . improved by M. Gurwell and B. But- 
ler (private communication, February 2002). R. Perley 
(private communication, February 2002) has shown that 
the flux density scale at 22.4 GHz based on this Mars 
temperature model agrees to within 0.4% with the flux 
density scale determined at this frequency based on a 
model of NGC7027 derived by A. van Hoof (private com- 
munication, February 2002). The ratio of the tempera- 
ture of Mars from the CBI to that of the improved Rudy 
model is 1.026 ± 0.013. Our second comparison is based 
on observations of 3C273 on the CBI and the VLA at 7 
epochs for which the time interval between observations 

2 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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on the two instruments was small enough to eliminate 
uncertainties due to source variability. The VLA obser- 
vations have been calibrated using the improved Rudy 
Mars model at 22.4 GHz, and the CBI 26-36 GHz mea- 
surements have been extrapolated to this frequency. The 
ratio of the flux density of 3C273 on the CBI to that on 
the VLA is 0.984 ±0.010. From the weighted averages of 
these two te sts we find that the CBI flux density scale, 
based on the lMason et alJ l|1999|) absolute calibration, is 
0.1%±0.8% lower than the VLA flux density scale at 22.4 
GHz. The level of agreement between these two scales is 
clearly fortuitous, but it does give us confidence that the 
5% calibration uncertainty we estimate for the CBI flux 
density scale is, if anything, conservative. 

Individual scans of calibrators and the microwave back- 
ground are bracketed by measurements of an internal 
noise calibration source whose equivalent flux density for 
each baseline and channel is determined by reference to 
the celestial flux density calibrators. The noise calibra- 
tion source was initially intended to remove instrumental 
gain fluctuations over the course of a night, but we found 
that the instrumental gains were more stable (~ 3% rms 
variations) than the noise calibration source output and 
so this was not used. This is primarily because of the 
difficulty in stabilizing the temperature of the noise cal- 
ibration source amplifier and various other components, 
all of which lie outside of the cryogenic dewars. All of 
the (complex) noise calibration source measurements for 
a given night are averaged together and the data from all 
baselines are scaled to give an identical response. This re- 
moves baseline-based gain and phase calibration errors; it 
also introduces antenna-based errors which are removed 
by the subsequent (antenna-based) primary flux density 
calibration. At the start and end of each night of observ- 
ing the relative gain and phase errors between the real 
and imaginary channels were measured using the noise 
calibration source with and without a 90° phase shift ap- 
plied to the receiver local oscillator. The rms quadrature 
phase is ~ 5° and the rms gain errors are ~ 10%. The 
solutions are stable over a timescale of several weeks. 

Lead and trail scans of the CMB fields are interleaved 
with 1 minute observations of bright (> 1 Jy) nearby cal- 
ibrators which provide a check on the telescope pointing. 
From these observations we determine that the absolute 
rms radio pointing is ~ 22", while the rms tracking errors 
are - 2". 

Data are flagged by both automatic and manual fil- 
ters. The on-line control system identifies corrupt or 
potentially unreliable data. These affect periods when 
the telescope was not tracking properly, a receiver was 
warm, a local oscillator was not phase locked, the total 
power of a receiver was out of the normal range, or when 
a receiver phase-shifter was not acquired. Observer notes 
were also used to examine periods when there were in- 
strumental problems that might have affected the data 
quality or there was bad weather, identified by visible 
cloud cover or corrupted visibilites on the short base- 
lines. A small fraction (~ 2%) of the data were deleted 
manually on the basis of these inspections. Subsequent 
automatic data edits are based on the data scatter and 
baseline-to-baseline correlations. The first level of rejec- 
tion is a 5cr scan-wise outlier edit to eliminate occasional 
instrumental glitches which give rise to large signals in 
individual samples. This cut affects a negligible fraction 



of the data, and our results are not sensitive to the pre- 
cise level of this cut. 

Further automated data-filtering is provided by our dif- 
ferencing procedure, which also generates an accurate es- 
timate of the thermal noise in each scan for each baseline 
and frequency channel. To do the differencing, individ- 
ual 8.4-s integrations in the lead and trail scans, which 
are taken 8 minutes apart, are matched and subtracted. 
Since our fields are separated by 8 min in observing time, 
the typical scan length is ~ 400 s after slewing and firing 
the noise calibration source, corresponding to ~ 50 inte- 
grations of 8.4 s each. In order to prevent a few short 
scans from biasing the statistics, any scan with fewer 
than 30 pairs of 8.4 s integrations was rejected. The 
noise estimates are derived from the differences between 
the lead and trail fields for individual scans of ~ 50 
8.4 sec integrations. Data from different baselines and 
times which contribute to the same (it, v) points are sub- 
sequently combined. This leads to an underestimate of 
~ 8% in the variance (see Appendix) which we correct for 
in our statistical analyses and spectrum determinations. 

The celestial signals are stationary on an 8-minutc 
timescale, and constant ground signal is cancelled by 
the lead — trail differencing. The noise will be increased 
by atmospheric noise, ground signal variations, or pre- 
viously unflagged instrumental anomalies. On this basis 
any scan for which the rms of the samples is greater than 
twice the expected thermal noise level is rejected: this 
primarily eliminates periods affected by weather. Less 
than 0.1% of the data are rejected by this filter. The 
noise level is established with reference to data which are 
clearly free of atmospheric contamination; it is highly re- 
peatable and within ~ 20% of the expected thermal noise 
level. 

In computing the data differences and thermal noise 
level, we also compute the covariance between all 
baseline-channel pairs, i.e., the thermal noise covari- 
ance matrix. For the deep 08h field, a single night 
(2000 January 12) was found to have statistically sig- 
nificant off-diagonal correlations. On this night there 
were clear instances of atmosphere emission on the short 
baselines, and the baseline correlations are likely due to 
poor weather. This day was excluded from the analy- 
sis. No anomalous correlations are seen in the deep 14h 
data. The correlation test on the deep 20h data revealed 
a single baseline-channel with a large real-imaginary anti- 
correlation. When this baseline-channel was eliminated 
from the 20h dataset, no significant baseline-baseline cor- 
relations are seen. The excision of these data from the 
08h and 20h fields has a negligible effect on the measured 
power spectrum, suggesting that lower-level, undetected 
correlator problems of a similar nature also do not affect 
our result. 

The telescope primary beam has been measured accu- 
rately using observations of Tau A. The primary beam 
width (FWHM) is 45.2 x (31 GHz/u) arc minutes; we es- 
timate a 0.7% uncertainty in the beamwidth. This gives 
rise to a 1.4% contribution to our overall calibration er- 
ror. Note that unlike the case for total power experi- 
ments, uncertainties in the primary beam width for an 
interferometer affect only the overall temperature scale, 
and do not introduce an ^-dependent systematic error. 
The beam- width measurement is described in Paper III. 
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Fig. 2. — CBI joint Deep Field window functions. The alter- 
nate binning window functions are displayed in the upper panel; 
the primary binning window functions are displayed in the lower 
panel. Vertical lines denote the bandpower bin divisions discussed 
in the text. The window function of the high-^! rebinning is shown 
as a dashed line. The expected value of band power for a given 
spectrum, C;, in one bin is E; W(l)/l X 1(1 + l)Ci/2ir 



4.2. Maximum Likelihood Analysis 

4.2.1. Formalism and Algorithms 

The determination of the angular power spectrum of 
the sky fr om interferometer me asur ements has b e en dis - 
cussed bv lHobson et all l|1995j) and lWhite et alJ ()1999|) : 
the more general fr amework of power spect rum estima- 
tion is discussed by iBond et al.l l)1998[) and iBond et al.l 
( 2000) . Here we present a brief overview of our proce- 
dure. The details of the method are given in Paper IV. 
We test a space of hypothetical models for the spectrum, 
Ci, against the data using the likelihood function as a fig- 
ure of merit. For complex data points V with zero mean 
and Gaussian noise described by the covariance matrix 
Cjk, this is 

L ^ = ^telC «P{-™-*V}. (1) 
where the covariance matrix, C, is given by 

C = C N + ]T q B C% + q src C src + q rcs C rcs . (2) 
B 

Here C N is the noise covariance and is the sky 
(or theoretical CMB) covariance for band B. The fore- 
ground covariance is made up of C src , a constraint ma- 
trix for point sources with known positions, and C res , 
a matrix which models the contribution of faint point 
sources of unknown position. The calculation of the ex- 
pecte d sky variance C R is discussed in many references 
(e.g-. lHobson et ai1ll995t Paper I; Paper IV). For a well- 
designed interferometer the noise term C N is diagonal 
(see also § I4.1f) . The pre- factors q rcs and q src are held 
fixed at a priori values (see S I4.3.2f> . As discussed below, 
the qs are varied to yield estimates of the microwave 
background bandpower. 

For the present analysis of the deep field data we adopt 
a parametric description of the Ci as seven bins with 
lower boundaries at I = 0, 500, 880, 1445, 2010, 2388, 
and 3000. Within each bin, the power spectrum is as- 
sumed to be flat in £(£ + l)Ce. The overall amplitudes of 



the bins qs are treated as free parameters in the maxi- 
mization of the likelihood. The spacing of A£ is set by 
the intrinsic resolution of the CBI in the aperture plane 
(A£ ~ 300). Here A£ is the FWHM of the visibility 
window function (see Paper III). This resolution can be 
increased by mosaicing. The divisions above were chosen 
to coincide with breaks in the multi-configuration (u, v) 
coverage, which minimize correlations between adjacent 
bins. Given the (it, v)-coverage of our data, this choice of 
binning yields 10-15% anti-correlations between adjacent 
bins. To demonstrate that our results are not sensitive 
to the choice of bins, we have also conducted the anal- 
ysis with an alternate set of bins whose boundaries are 
halfway between those of the primary binning described 
above (six bins with lower boundaries at £ = 0, 690, 1162, 
1728, 2199, and 2694). Window functions for each band 
are shown in Figure [21 These show the sensitivity of a 
given bandpower bin to power in individual multipoles 
as a function of £. 

We compute the parameter uncertainties and correla- 
tions from the curvature of the likelihood at the best-fit 
locus. In practice we use the Fisher matrix 

as a computationally efficient estimator of the curva- 
ture. In the approximation that the likelihood function 
is Gaussian, Fg B , is the correlation matrix of the fitted 
parameters. As a check on the Fisher matrix results, we 
have also mapped the likelihoods in individual bands and 
checked against the err or bars computed i n the offset log- 
normal approximation ( Bon d et alJl200Cjl . We find good 
agreement between these three methods. Monte Carlo 
tests of simulated CBI data show that these error esti- 
mates, as well as the bandpower estimates themselves, 
are unbiased. 

The dominant foreground for the CBI is composed of 
discrete extragalactic radio sources described by the cor- 
relation matrices C rcs and C SIC . For purposes of this 
analysis we neglect other possible foregrounds, although 
see § |5] for limits on the possible contribution of fore- 
grounds on large angular scales. 

4.2.2. Implementation 

Since in the small-angle limit interferometers directly 
measure linear combinations of the sky Fourier modes, 
the extraction of angular power spectra from the data is 
in principle simpler than in the case of total power exper- 
iments. In general, however, Cjk has non-negligible off- 
diagonal elements, and depending on the (it, v)-sampling 
of the dataset this matrix can be large. To deal with 
the computational challenge of inverting this matrix, we 
have developed a Fourier gridding algorithm. This algo- 
rithm constructs a data vector composed of linear combi- 
nations of the visibility data and calculates the resulting 
covariance matrix. For single-field analyses such as we 
report here, the primary effect of this gridding is to re- 
duce redundancies in the (u, w)-dataset caused by dense 
aperture plane sampling. This compression is typically a 
factor of 3 to 4: for each field there are several thousand 
visibilities, which are gridded down to 970 estimators. 
This results in over an order of magnitude increase in 
the speed of the likelihood analysis. This pipeline and 
the associated formalism are discussed in Paper IV. 
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The pipeline was implemented on the 4 and 32 pro- 
cessor shared memory Alpha GS320 and ES45 clusters 
at CITA. The codes were parallelized using Open MP 
directives. A complete joint analysis of the three deep 
fields takes approximately 4 CPU hours. This includes 
pre-gridding and the calculation of all ancillary quanti- 
ties such as the bandpower window functions, likelihood 
maps and noise spectra. Bandpower estimates can be 
obtained in as little as 1 CPU hour for the deep field 
data. These rapid turnaround times permitted extensive 
testing of the data and pipeline; results of some of these 
tests are discussed in i)4.4l 

4.3. Discrete Sources 

Extragalactic radio sources are the dominant fore- 
ground to CMB observations ov er the 26-36 GHz fre- 
quen cy band on arcminute scales ijTegmark fc Efstathioul 
1996). If no allowance for this foreground is made, the 
few brightest discrete sources in each field dominate any 
power spectrum determination beyond £ ~ 1000. While 
much is known about radio source populations at lower 
frequencies, this is not the case at 31 GHz, and observa- 
tions at this frequency are essential in dealing with the 
point source contamination at high £. We therefore con- 
structed a 26-34 GHz receiver for the OVRO 40 meter 
telescope which we have used to study the source pop- 
ulation at this frequency and to provide a direct check 
on the CBI calibration. We have also computed 31 GHz 
source counts directly from the CBI deep maps, and the 
mosaic maps presented in Paper III. The results from 
these complementary approaches are presented in § 14.3.11 
These results form the basis of our strategy for dealing 
with both the bright sources and the statistical residual 
source background. The treatment of radio sources in 
the CMB analysis is detailed in § 14.3.21 

4.3.1. 31 GHz Radio Source Measurements 

The OVRO 40-meter telescope was used at 31 GHz to 
survey the 2225 NVSS sources brighter than Si.4GHz = 
6mJy in four 22.5 square degree fields which encompass 
the deep and mosaic fields. With the 4<r cutoff speci- 
fied below, the survey has a 7% chance of one or more 
false detections. To minimize the effects of variations in 
the flux density of the monitored sources, the 40-meter 
and CBI observations were made as nearly as possible 
over the same period of time. The typical sensitivity 
achieved on the 40-meter telescope was 2 mJy (rms) , al- 
though the distribution extends to 5 mJy. 12 % of the 
targeted NVSS sources were detected above the survey 
thresholds of 4cr and S m in,3i = 6 mJy. The survey is 
90% complete at S 3 i > 16mJy and 99% complete at 
5*31 > 21mJy. The 56 detected sources within 1 degree 
of the lead or trail field centers were subtracted directly 
from the CBI data; these sources are shown in Table |5J 
This table lists total flux densities, S31; distances in arc 
minutes from the center of the CBI field, r; and the flux 
density predicted on the CBI including the effects of the 
CBI primary beam, Sqbi- F ew of these 56 sources con- 
tribute significantly to the data since most of them are 
too distant from the CBI field centers to do so: in each 
deep field only a few sources are visible in the images. 
The spectral indices of these sources were computed by 
comparison with the NVSS 1.4 GHz flux densities, and 
we find a — —0.45 with an rms dispersion of 0.37; the 



minimum and maximum spectral indices were —1.32 and 
+0.50. 

We have used the CBI deep and mosaic maps to de- 
termine the source counts at 31 GHz. To do this, we 
made maps from the long baseline (> 250 A) data and 
searched for peaks over the 5cr threshold. Below this 
level the probability of false detections is significant. A 
fit to the resulting counts over the range 5-50 mJy yields 

x -1.0 



A(> S31) = 2.8±0.7deg" 



•5': 



31 



(4) 



10 mJy / 

This is slightly higher than the 31 GHz counts estimated 
from either the NV SS counts or the 15 GHz counts of 
iTavlor et all (|2001^ . In both the deep and mosaic maps 
no sources were detected above the 5a limit which were 
not NVSS objects. The area searched depends on the 
limiting flux density chosen: in the deep maps the area 
was 6 square degrees at 12 mJy, falling to 1.8 square 
degrees at 6 mJy. In the mosaic maps the area was 47 
square degrees at 25 mJy, falling to 12.2 square degrees 
at 18 mJy. The flux densities of the sources detected 
in both the deep and mosaic maps are consistent with 
the OVRO determinations, with the agreement typically 
within 5%. 

As discussed above, the OVRO-detected bright sources 
were directly subtracted from the CBI data since the 
OVRO source subtraction is a step in the standard data- 
analysis pipeline which has been useful in producing diag- 
nostic images of the data. In § l4.3.2lwe discuss a method 
(the constraint matrix approach of H alverson et al.l2002j) 
of dealing with point source foregrounds which does not 
require knowledge of the source flux densities. We find 
that this method applied to the OVRO sources gives very 
similar results as direct subtraction. These complemen- 
tary techniques of dealing with this bright source popu- 
lation have been useful in extending our power spectrum 
measurments to the instrumental limit of 3500. 

4.3.2. Statistical Treatment of Sources 

In addition to the brightest sources in the CBI fields, 
which have been measured and subtracted with the 
OVRO data, we must deal with sources too faint to mea- 
sure directly but which contribute statistically to our 
measurements. Since most radio sources have spectra 
which fall towards higher frequencies, and since there 
are wide-area radio sky surveys at lower frequencies with 
sensitivities comparable to what we have in the CBI deep 
field maps, most of the sources which contribute are read- 
ily identified from the low-frequency surveys. For this 
purpose we use the NVSS. There is also a small contri- 
bution due to sources which are too faint to appear even 
in NVSS. We call these latter "residual" sources and es- 
timate their contribution from the CBI 31 GHz source 
counts. 

After the subtraction of the bright sources monitored 
at OVRO, there remain three components of point source 
contamination: one due to data correlations induced by 
the (imperfect) OVRO source subtraction; a component 
due to the NVSS sources whose positions are known, but 
whose flux densities are not measured at 31 GHz; and 
a residual component due to sources whose contribution 
must be estimated from source counts. These are repre- 
sented as 

pOVRO , „ n NVSS , „ ^res ft r\ 
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Table 2. Subtracted Discrete Sources 
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The q are prefactors which are held fixed in the band- 
power fit, and allow for the effect of overall uncertainties 
in our estimates of the source covariances. If q = 1, an 
extra amount of variance equal to the variance assumed 
in the calculations of the source covariance matrices is al- 
lowed in each source mode. In the limit q — > oo the data 
modes corresponding to the point sources are completely 
remo ved from the data , i.e., these modes are "projected 
out" l|Bond et al.lll998HHalverson et alJ l2002). We note 
that it is possible to allow the q to be free parameters in 
the fit, and we do not find significantly different results in 



this case; however it is more conservative to fix the q at 
very large values. Explicit procedures for calculating the 
three source covariances matrices are discussed in Paper 
IV. 

We adopt the strategy of projecting the modes corre- 
sponding to OVRO residuals and NVSS sources out of 
our data. A nominal variance is assigned to each source 
based on the OVRO measurement uncertainty or, for 
the NVSS sources, the width of the OVRO-determined 
1.4 to 31 GHz spectral index distribution. We then set 
Qovro = Qnvss = 30 which projects these modes out of 
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the data. We find that our power spectra at high £, where 
point sources are most troublesome, are little affected by 
variations in q between q = 30 to q = 10 6 . 

At significantly larger values of q the matrices are nu- 
merically ill-conditioned. We have tested the source pro- 
jection algorithm extensively on both real and simulated 
data and find that it is robust against, for instance, ran- 
domly reassigning estimated variances between sources 
and randomly perturbing the source positions assumed 
in the projection by 30" (rms). 

We model the contribution of the residual sources 
(C res ) as a Gaussian white-noise foreground. We de- 
termine the normalization of this matrix from Monte 
Carlo simulations assuming the source counts presented 
in § 14.3.11 to generate random populations at 31 GHz. 
Since at fainter flux levels the counts are expected to 
flatten somewhat, we used a slightly shallower slope for 
the counts (N(> S) oc S' - 875 ) for this calculation, 
although at the la level our result is consistent with 
that obtained for N(> S) oc S^ 1 . In order to deter- 
mine what fraction of these sources fall below the NVSS 
detection threshold of 3.4 mJy and hence are residual, 
the sources are extrapolated down to 1.4 GHz using the 
observed OVRO-NVSS distribution of spectral indices 
(with a = —0.45 — see § 14. 3. II) . Since the OVRO survey 
is strongly biased towards flat-spectrum sources this will 
overestimate the residual source power level. We find 
a 31 GHz variance of Cl cs = 0.08 Jy^r" 1 . We get a 
similar result, Cg es = 0.10 Jy 2 sr _1 from the analytic cal- 
culation described in Paper IV. In our highest bin this 
corresponds to a power level of ~ 115 /iK 2 , which is less 
than the thermal noise in this bin and all others; we as- 
sign an uncertainty of 50% to this value, corresponding 
to ±0.04 Jy 2 sr _1 (57 /iK 2 ). This uncertainty is based on 
varying the number counts and spectral index parame- 
ters within the range allowed by their uncertainties, and 
rerunning the Monte Carlo simulations used to calcu- 
late the residual power. Comparable power levels are 
obtained if, instead of using the 31 GHz counts, we ex- 
trapolate in frequency from the 1.4 GHz number counts 
using a simple Gaussian model for the spectral index dis- 
tribution which is constrained to reproduce the detection 
statistics observed in the OVRO survey. As a test on the 
robustness of this correction we re-ran the power spec- 
trum analysis doubling the NVSS/rcsidual source thresh- 
old from 3.4 to 6.8 mJy at 1.4 GHz with similar results. 
The number density of the residual sources on the sky is 
too high to permit full projection of these data modes, 
and in any case their positions are not known so this 
would be impossible. We therefore fix g ros = 1, which 
approximately corresponds to subtracting our conserva- 
tive estimate of the residual source power level from the 
spectrum. 

The cost of projecting out the NVSS sources is a loss of 
sensitivity of the CBI at high I due to our limited knowl- 
edge of the faint point sources at 31 GHz. The best way 
to address this problem in the near term is with a sen- 
sitive 31 GHz receiver on the NRAO Green Bank Tele- 
scope (GBT). Such an instrument is under construction 
and should be in use by the end of 2003. 

4.4. Data Consistency Tests 

The redundancy of observations on the CBI deep fields 
allows many checks on the dataset. Below we describe 



the key tests performed. 

4.4.1. x 2 Tests 

A simple and direct test of the data consistency is pro- 
vided by computing the x 2 °f visibilities with the same 
(u, v) coordinate, frequency, and sky pointing. We com- 
pared the data at each (u, v) point on each night with 
the average of all other nights. This procedure identified 
three days on the deep 20h field which were formally in- 
consistent with the total dataset. We exclude these days 
from the analysis, which has a small (<C la) effect on our 
results. The visibilities for all other days were consistent 
with the estimate of the visibilities from the rest of the 
data. A typical value of the x 2 °f an individual day as 
compared with the rest of the dataset is x 2 l v — 0.986 
with v = 1617. For such a day the probability of exceed- 
ing this value under the null hypothesis is 65%. 

The data on each configuration of each field was also 
subdivided in half by time, subtracted, and then com- 
pared in this manner. The x 2 l v values for these tests 
on the 08h data (3 configurations), and the 14h, and 20h 
data (each one configuration) were 0.990, 0.993, 0.967, 
1.015, and 0.986 for v = 2310, 2710, 3948, 2310, and 
4848, respectively. These results indicate that our esti- 
mate of the thermal noise variance is correct to within 
~2%. 

4.4.2. Kolmogorov-Smirnov Tests 

We have used the Kolmogorov-Smirnov (KS) test to 
determine if the visibility data are consistent with a 
Gaussian distribution. This is the expectation when the 
signal-to-noise-ratio per baseline-channel combination is 
low, as is the case for the I > 1000 data on any given 
field once the bright point sources have been subtracted. 
The quantity of interest is the ratio of the visibility to 
its estimated uncertainty 



If the ay were known a priori then t would have a Gaus- 
sian distribution with unit dispersion. As discussed in 
the Appendix, however, this is not the case, and a 4% 
correction to the expected standard devia tion (a t ) of t is 
required. Furthermore, as discussed in § 14.4.11 we have 
only determined that our measurements of the thermal 
noise are accurate at the 1% level. We therefore apply 
the KS test to the data under the null hypothesis that 
the observed distribution of t is consistent with a Gaus- 
sian distribution of zero mean and at — 1.04, but vary 
at by ±1% for each dataset. We find good consistency 
with the Gaussian distribution, so to the accuracy that 
the noise level has been well-quantified we do not detect 
significant non-Gaussianity in the visibility data. 

4.4.3. Power Spectra 

Two sets of bandpowers q\ t g an d q2,B (i-e., two power 
spectrum estimates) can be compared for consistency by 
forming 

X 2 = (qi,B> - <72,B') T (FiBB> + F 2,BB')~ l («1,B - Q2,b)- 

(7) 

Under the hypothesis that the likelihood function is 
Gaussian and the qi.B are drawn from the same underly- 
ing power spectrum with noise correlations specified by 
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the elements of the inverse of the Fisher matrices F i BB ,, 
this is distributed as \ 2 with v degrees of freedom, where 
v is the number of free parameters in each spectrum es- 
timate (7). The x 2 l v values (and significances) for the 
08h/14h, 08h/20h, and 14h/20h power spectrum com- 
parison are, respectively, 1.40 (20%), 1.03 (40%), 1.08 
(38%). 

4.4.4. Undifferenced and Doubly-Differenced Power 

Spectra 

In order to estimate the level of non-celestial signals 
which we are removing from the data, we have extracted 
the power spectra of the undifferenced data (i.e., no 
lead — trail based ground subtraction) . The power spec- 
trum estimates derived from the shortest baselines show 
the effect of ground-spillover, but at I > 600 they are 
within a factor of ~ 4 (in /iK) of the differenced spec- 
trum, and are seen to fall strongly with increasing £. 
Since the CMB signals on the short baselines are up to 
a factor of ~ 10 3 (in /iK) less than the ground signal on 
these baselines, the low power levels in the undifferenced 
spectrum indicate that the signal averages down when 
data from many days and position angles are combined. 
At multipoles £ > 1000 the undifferenced signal is within 
~ ler of the differenced spectrum except in a single bin. 
In this bin (2388 < i < 3000) the undifferenced power 
level was originally found to be about twice the thermal 
noise level. This was traced to correlator offsets in a 
single channel on several baselines. Since the x 2 of the 
differenced data on these baselines are consistent within 
the noise with data from other baselines and the power 
spectrum is not significantly affected by the exclusion of 
the channel in question, we infer that the contamination 
is effectively removed by our differencing procedure. 

We have also divided the data from the lead and trail 
pointings for each field in half, subtracted these, and 
derived the power spectrum of the combination of the 
doubly-differenced data sets for the three fields. This is 
consistent with zero — X 2 l v — 1-06 for v = 7 DOF (7 
bins). 

4.4.5. Simulated Datasets 

We carried out a series of Monte Carlo tests on the 
data. In these tests, the individual days that went into 
the deep field datasets for which we have high sensitiv- 
ity, 08h and 20h, were resampled with the observed (u, v) 
coverage and noise assuming a known input power spec- 
trum. Since the signal to noise ratio beyond i ~ 1000 in 
a single day is very low, we extracted only the first two 
bandpowers from each day and its corresponding simu- 
lated dataset. We find that the distribution of the band- 
powers of the real data is consistent with that of the 
simulated data. The distribution of the rms powers of 
the real data within the primary beam area on the dirty 
maps is likewise consistent with that of the simulated 
data. 

5. RESULTS 

The power spectra of the individual deep fields are 
shown in Figure 01 These illustrate the field-to-field con- 
sistency of power spectra discussed in § 14.4.31 Figure 01 
shows the results for the joint fields, as well as the ther- 
mal noise and residual source power spectra. It is evident 




Fig. 3. — Power spectra from the individual CBI deep fields 
shown for the primary binning only. The 08h, f4h, and 20h data 
are denoted by squares (blue), circles (green), and crosses (red), 
respectively. The horizontal location of each point on the power 
spectrum is £ e fj and the horizontal error bars indicate the full ex- 
tent of each band. 



He 




1000 2000 3000 4000 

Fig. 4. — Joint CBI deep field power spectra for the primary 
(blue open circles) and alternate (green open squares) binnings. 
The placement of the points in I is determined as in figure Also 
shown are the thermal noise power spectrum (open stars) and the 
residual source power spectrum (solid triangles). 



that the residual source correction is small compared to 
the observed power levels, and that errors in the thermal 
noise level power spectrum at the 1% level indicated in 
§ 14.4.11 are also small compared to observed power lev- 
els over the entire ^-range. The bandpower results are 
summarized in Table |3 In this Table, £ e ff denotes the 
centroid of the bandpower window function. The error 
bars in this table are la 

Figure compares the CBI deep field results directly 
with the BOOMERANG, DASI, and MAXIMA data. 
We see that there is good agreement between all the 
observations in the range i ~ 600 — > 1000, but that 
the CBI deep field observations are somewhat higher 
than the other experiments in the range £ < 600. This 
discrepancy — which is not statistically significant — is 
discussed further below and in § We see in this fig- 
ure that the power level drops significantly in the range 
£ ~ 500 — » 2000, confirming the results of Paper I, which 
were based on a subset of these data, and extending the 
region in ^-space over which this drop occurs to higher 
£. Above £ ~ 2000 the level of power is flat and signif- 
icantly greater than zero, contrary to what is expected 
from standard models of the intrinsic anisotropy. We 
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Table 3. Deep Field Bandpower Data (fiK 2 ) 



Bin 






Joint 


deep 08h 


deep 14h 


deep 20h 


- 500 


307 


7193 ± 2478 


6531 ±4100 


8381 ± 5274 


6829 ± 3540 


500 - 


- 880 


640 


1896 ± 556 


3213 ± 1306 


910 ± 983 


1322 ± 627 


880 - 


- 1445 


1133 


724 ± 200 


1045 ± 389 


310 ±461 


674 ± 264 


1445 


- 2010 


1703 


245 ± 149 


449 ± 266 


-415 ±516 


578 ± 257 


2010 


- 2388 


2183 


657 ± 323 


251 ± 382 


1174 ± 1620 


1293 ± 557 


2388 


- 3000 


2630 


573 ± 308 


795 ± 442 


2742 ± 1897 


-5 ± 320 


3000 


- 4000 


3266 


-37 ±418 


-276 ± 435 


3717 ± 3195 


604 ± 920 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 - 




S^MS BOOMERANG 
DASI 
MAXIMA 



-e- 



1000 



2000 



3000 



4000 



I 



Fig. 5. — Comparison of the C BI deep field power spectrum with BOOMERANG iNctterficld ct al. 200$, DASI I Halvcrso n et all2002fl . 
and MAXIMA ILee et al J 1200 II) results. The rectangles indicate the 68% confidence intervals on band-power; for BOOMERANG, the 
solid rectangles indicate the 68% confidence interval for the statistical and sample variance errors, while the hatched rectangles shows the 
amount by which a ±1<t error in the beamwidth (12.'9± 1.'4) would shift the estimates (all up or all down together). Results from the CBI 
primary binning are shown as blue open circles, and the alternate binning results are shown as green open squares. The high value of the 
power seen in the first bin of the primary and alternate binnings relative to other observations has rather low significance (~ 1.7c) due to 
the large sample variance in the deep field spectra at low i\ this is discussed further in the text. 



return to this point in §5.2. 

In Figure |HJ we have plotted the CBI deep spectrum 
together with the best-fit model to the CBI mosaic re- 
sults combined with the BOOMERANG, DASI, MAX- 
IMA, VSA and earlier observations ( "all-data" ) , derived 
in Paper V, subject to the constraints listed in the figure 
caption. Here we also show the expected level of the CBI 
deep spectrum in the primary and alternate binnings ob- 
tained by integrating this model over the CBI window 
functions shown in Figure |21 We find that the CBI deep 
spectrum is reasonably consistent with the theoretical 
curve at £ < 1000 — the first two bins of the standard 
binning yield \ 2 jv = 1.86 (84% significance) to the the- 
oretical curve, and the first two bins of the alternate 
binning yield x 2 /^ — 2-03 (87% significance). 

We have constructed Wiener-filtered images of the sky 
signal in which we have subtracted the point source and 
noise contributions by th e method described in Paper IV 
andin lBond et all lfl9941. The images for the three fields 
are shown in Figured 



5.1. Limits on Diffuse Foregrounds 

We can use the spectral coverage of the CBI to place 
limits on potentially contaminating foregrounds such as 
Galactic synchrotron (a ~ —0.7) and free- free emission 
(a ~ 0). To do this, we generated a set of 100 realiza- 
tions of our data set using the best-fit power spectrum 
shown in Figure H3 as well as 100 realizations of a fore- 
ground with zero spectral index, typical of free- free emis- 
sion, with a flat power spectrum. We added the two sets 
together at various power levels of the foreground, and 
fitted the data from I < 800, where we are most sensi- 
tive to the spectral index, to a two-bin ( £ < 500 and 
£ > 500) model with a varying spectral index, finding a 
single best-fit spectral index using that model. Models 
containing both CMB and foreground components are 
too ill-constrained by the data to be useful. We find that 
the mean best-fit spectral index for the zero-foreground 
simulations is 1.97, and that the scatter of the individual 
simulations is 0.34. Furthermore, we find that the mean 
best-fit spectral index at different simulation power lev- 
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* Model bandpowers 

CBI high I band (68%, 95%) 




1000 



2000 
I 



3000 



4000 



Fig. 6. — The CBI deep field power spectrum: Results from the CBI primary binning are shown as circles, and the alternate binning results 
are shown as squares. The shaded region shows the 68% and 95% central confidence intervals for the high-^ CBI bin 2000 < £ < 4000; these 
was computed by combining the last three primary bins into a single bin in bandpower and mapping the likelihood of this band, marginalizing 
over irrelevant parameters. The curve shows the spectrum derived in Paper V as a best fit to "all-data" (Boomerang, CBI, DASI, DMR, 
Maxima, VSA and earlier observations) using strong priors (age of the universe > 10 Gyr, 45 kms -1 Mpc -1 < Ho < 90 kms" 1 Mpc - 1 , 
zero curvature, plus large scale structure constraints). Stars show the expected signal in each CBI band for this model. The power in the 
range I ~ 2000 — > 3500 differs from the best fit model at the 3.1<r level. This excess power is discussed further in the text. 
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mjy/beam 
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Fig. 7. — Wiener-filtered images of the 08h, 14h and 20h deep fields (see § 8 of Paper IV). The positions of the field centers are given in 
Table 1. The dominant features in these images have angular scales ~ 15', and are due to structures with multipoles observed on ~ 1 meter 
baselines, which span the range of the second and third acoustic peaks. The approximate fraction of the data contributing to these features 
can be seen in Figure^ where the peaks in the range 500 < I < 800, correspond to the data contributing to the dominant features in these 
images. Due to the differencing and the sidelobes of the point spread function, there is some ambiguity in relating particular features in 
these raw maps to actual features on the sky. Figure 1 of Paper I presents a similar image of the 08h field, along with the point spread 
function and one of several possible deconvolutions. 



els is well described (1% error) by taking the means of 
the CMB and foreground spectral indices, weighted by 
their power levels. 

^cmb-^cinb C^fgPfg , * 

"fit = p T-75 . (») 

r cmb "T -* fg 

with P cmb = 2455 ^K 2 and a cmb = 1.97 for 500 < t < 
800. The data give a best-fit spectral index of 1.84±0.34, 
0.4cr from the zero-foreground value. The la upper limit 
for I < 880 on a free-free like foreground contribution 
is aet = 1.84 - 0.34 = 1.50, which, using EqQ3 yields 
P fg = 744 fiK 2 (31% of P cmb ) for a fg = 0. For a syn- 



chrotron spectral index of —0.7, we get an upper limit 
Pfg < 521 nK 2 (P fg < 21% of P cm b). 

5.2. The Apparent Excess Power at High £ 

An interesting feature of the joint spectrum is the ap- 
parent excess power observed above I — 2000. To quan- 
tify this, we reanalyzed the data with the last three 
bins in the primary binning grouped into a single bin 
(2010 < t < 4000). For the single high-^ bin we ob- 
tain a best-fit bandpower of 508 /zK 2 with an uncer- 
tainty (from the Fisher matrix) of 168 /iK 2 . In order to 
more accurately estimate the uncertainty we have calcu- 
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Table 4. Effect of Source Treatment on High-<? Bin 



Source Corrections 



Power Level 
QK 2 ) a 



H850 
-718 
M4!> 
-361 
j+372 
-279 

OVRO, NVSS, and Residuals 508 ± 356 



No Correction 

OVRO Subtraction Only 

OVRO+NVSS Corrected 



3028Z 
1074^ 
603 4 



a Errors are ±2<r. 



lated the likelihood curve of the bandpower in the high- 
t bin explicitly — marginalizing over the bandpower in 
the previous bin, which is the dominant bandpower cor- 
relation, and the residual source correction — and inte- 
grated the likelihood function to determine confidence 
intervals. The 68% and 95% central confidence intervals 
are, respectively, (359 -> 624) /zK 2 and (199 -> 946) /zK 2 . 
Averaged over the high-^ bin, the model shown in Fig- 
ure predicts ~ 66 /iK 2 of power. Power levels this low 
or lower are excluded at a significance corresponding to 
3.1<7, and zero signal is excluded at a significance corre- 
sponding to 3.5(7. The best-fit values of the other four 
bins of the high-£ analysis are close to those achieved 
in the primary (seven-bin) analysis. The joint deep-field 
power spectrum is shown with the high-£ 95% central 
confidence interval (~ ±2er) in Figure [5] These results 
differ from what would be obtained from the Fisher ma- 
trix error bars since — although the Fisher matrix calcu- 
lates the overall curvature accurately — the actual like- 
lihood is asymmetric. Using the Fisher matrix uncer- 
tainties the significance of a detection of non-zero power 
at I > 2010 is 2.8cr, including the contribution of the 
uncertainty in the residual source power. 

In order to give a clear idea of the corrections being 
applied, Table 0] shows the changes in the high-£ bin 
that result from the application of the individual cor- 
rections. Two-thirds of the source power is eliminated 
by the OVRO subtraction. Of the remaining ~ 500 [iK 2 
of excess source power, about 80% is removed by project- 
ing NVSS sources out of the data, and 20% is removed by 
the statistical residual source correction. As indicated in 
S I4.3.2l thc known source projection is robust; the results 
are insensitive to even large changes in the projection 
coefficients, q, and to fairly major corruptions of the co- 
variance matrices such as randomly reassigning variances 
between sources. 

The spectral index distribution of the OVRO sources, 
used in calculating the residual source correction, cov- 
ers the range —1.3 < a < +0.5. However, sources 
with spectral indices up to a = +2 have been detected 
Ce.g. JEdge et alJl998l:IStanghellini et al.ll998|) . We have 
therefore explored the possibility that a seperate popula- 
tion of objects, not seen in the NVSS or accounted for in 
the residual source correction, might be responsible for 
the excess. For this hypothetical population, we assume 
a power-law integrated number-flux density slope (3 at 31 
GHz, and we then compare the number density of sources 
in this hypothetical population required to produce the 
observed excess with the limits from the CBI data, where 
we found no sources above the 5-cr cutoff in the deep 
or mosaic fields which were not correctly identified with 



NVSS sources, as discussed in § 14.3.11 Even for the very 
steep integrated source count slope of /3 — —1.5, which 
has not been observed at low flux densities and high fre- 
quencies, we find that the number density of sources in 
this hypothetical population required to explain the ex- 
cess exceeds the upper limit from the CBI observations 
by > 3cr. It is therefore unlikely that such a hypothetical 
population of inverted spectrum sources is responsible 
for the excess. 

The apparent excess would be explained if we had un- 
derestimated the residual source correction by a factor 
of 4.5. We have been unable to construct a model which 
achieves this while remaining consistent with the source 
counts we have derived at 31 GHz, and source statistics 
at other frequencies. 

6. DISCUSSION AND CONCLUSIONS 

In this paper we have presented measurements of the 
CMB power spectrum out to I = 3500, beyond the scales 
probed by BOOMERANG, MAXIMA, and DASI and 
well into the damping tail region of the spectrum. Our 
present results confirm our previous detection of a drop 
in power at multipoles above I ~ 750 relative to the 
level at lower multipoles, based on a subset of the data 
used here (Paper I), and show that the decline in power 
persists out to I ~ 2000. Such damping is one of the 
fun damental p redictions of standard cosmological mod- 
els tfSiTlJ 11968). Below I ~ 900 the power levels ob- 
served are consistent with those seen in earlier exper- 
iments dde Bernardis et all |2000 t IHalvers on et all 12002 : 
Hananv et alJl2000t iLee et all 1200 H iLeitch et alJboOC : 
Miller et al.ll999fl . At I < 500 the power detected here is 
greater than that seen in other experiments including our 
own mosaic observations (Paper III). As discussed in §|SJ 
most of this discrepancy is due to the first bin. The CBI 
mosaic spectrum (Paper III), which is based on a larger 
area and therefore has lower sample variance than the 
deep field results, has power levels at low t which show 
excellent consistency with the BOOMERANG, DASI, 
and MAXIMA results. It is therefore likely that the dis- 
crepancy seen in the deep fields is due to the large sample 
variance in the \ow-£ mode estimates from these fields. 

In Paper V we discuss the constraints which the CBI 
deep and mosaic data provide on standard cosmologi- 
cal parameters. The CBI mosaic fields are much more 
powerful for this purpose than the deeps due to their 
higher ^-resolution and lower cosmic variance. The deep 
fields, however, provide a robust check on the results: 



we find fit 



o.6i: 



■0.10' 



Hot = l.Q9l8:&; n. = 1.16±g;i|; h 
and an age for the universe of 12.2 ±1.8 Gyr. These re- 
sults assume the weak-/i prior (discussed in Paper V) and 
use only the deep field power spectra out to I = 2000, 
where consistency with standard models of the instrinsic 
anisotropy is observed. More discussion of these results 
and the analysis method can be found in Paper V. 

Above I ~ 2000 we detect a 3.1c excess in power rel- 
ative to the best-fit curve. Possible explanations of this 
excess are: 

1. Data analysis error. We have tested our analysis 
by developing two independent software pipelines, 
and by developing a simulation program which gen- 
erates data sets which mimic the (u,v) coverage 
and the distribution of measurement uncertainties 
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of the real data set precisely. We have subjected 
both real and simulated data analyzed with both 
pipelines to a large battery of tests, and have been 
unable to find significant inconsistencies. 

2. Instrumental problems. Any instrumental signals 
which are stable over an 8-minute time span are 
removed to high precision by the differencing; in- 
strumental signals which vary on < 8 minute time 
spans would be evident in the doubly-differenced 
power spectra, but we see no such signals. It is 
also possible that there are inadequately modeled 
instrumental effects. A prime candidate would be 
pointing errors, but our pointing errors of ~ 20" 
are too small to account for an effect of this mag- 
nitude. 

3. Primary CMB anisotropy. This is inconsistent with 
standard theories which fit the low-£ range well, 
and is therefore an unlikely explanation. 

4. Diffuse Galactic foreground. We cannot rule 
out this possibility with the present observations. 
Higher sensitivity observations at the same fre- 
quency or at a higher frequency could test this 
possibility. We cannot account for the signal with 
known diffuse foregrounds, and the sensitivity is 
too low to constrain the spectral index of the sig- 
nal at these multipoles. However, in view of the 
anomalous component of Galactic e mission that 
has be en detected at this frequency bv lLeitch et al.1 

we are pursuing correlation analysis of the 
sky images with the IRAS 100 /im flux density, and 
with other signals, such as the Ha intensity. 

5. Residual point sources. This appeared initially 
to be a likely candidate, but the constraint ma- 
trix approach to removing point source foregrounds 
has proved to be remarkably robust, and our 31 
GHz source counts have enabled us to place strong 
limits on the hypothetical population of inverted- 
spectrum point sources needed to produce the ex- 
cess. While we cannot strictly rule out such a pop- 
ulation, a very steep integral counts slope would be 
implied (7 = —1.5 or steeper), as well as a normal- 
ization which is inconsistent at the 3cr level with 
that determined from the CBI-determined source 
counts. Future 30 GHz surveys with the GBT will 
allow this issue to be further addressed. 

6. Secondary anisotropy. There has been great 
interest in predicting the nature of the statis- 
tical SZE contribution to the CMB anisotropy 
on s mall angular scales, using bo th analytical 
(e.g., iCooravlhoOlt iMa fc FrvllgOOl and numer- 
ical (e.g.. IBond fc Mversl 119961 IBond et al] 12002]: 
Ida Silva et all 120001: iRefreeier fe Tevssierl booi 
ISeliak et afl T2002) methods. These works gener- 
ally predict a crossover between the intrinsic CMB 
and SZE signals at £ ~ 2000-3000. The level of 
SZE anisotropy forecast by theoretical models is in 
the range of 30 /1K 2 to 300 /iK 2 , depending sensi- 
tively on the rms mass fluctuation on large scales in 
the present universe (characterized by erg). There- 
fore secondary SZE anisotropy is, at some level, 



likely to contribute to the excess we report, but it 
is not clear if the majority of the observed signal 
can be attributed to the SZE. If SZE anisotropics 
were to be the cause of the observed excess, values 
of erg > 1 would be favored. For a detailed dis- 
cussion of the possible implications of the observed 
excess for models of SZE anisotropies, see Paper 
VI. 

Other possible contributors to signals on these angular 
scales include t he Vishniac effec t llVishniacll98^. patchy 
re-ion ization l)Aghanim et alJ 119961: iGruzinov fc Hul 
1998) and gravitational lensin g iBlancha rd fc Schneider! 



lyyoj and gravitational lensing (rilancnard &l bchneider 
T^ICoIefc EfstathiodlllMl^elirklll996j) . All of these 
are expected to be small effects compared to the signal we 
observe. It should be borne in mind that if the signal is 
due to non-Gaussian structures then the sample variance 
errors in our result and others will have bee n underesti- 
mated; see, for example, iZhang et alJ l)2002ft and Paper 
VL 

iDawson et alJ l)2001[) have reported a tentative (1.3cr) 
detection with the BIMA array at the same frequency 
on a smaller angular scale. These investigators place a 
95% upper limit of Qfl at = 12.4 ^K, corresponding to 
t{l + l)Q/2?r = 369 /iK 2 at £ cS = 5530. Since this mea- 
surement is at higher £ than the excess we have found 
here, it is possible that these two results are not directly 
comparable, but as discussed in Paper VI they could both 
be affected by secondary SZE anisotropy. 

The key result of this paper is the clear demonstration 
of the existence of a damping tail to the anisotropy spec- 
trum over the range £ ~ 500 — > 2000. This shows that 
on average there are no large deviations of the intrin- 
sic anisotropy spectrum from the predictions of standard 
cosmological models over this £-range. As discussed in 
Paper III, these measurements also support the gravi- 
tational instability paradigm for structure formation in 
the universe by providing the first direct measurements 
of the seeds from which present-day galaxy groups and 
clusters formed. In addition we report a detection of 
power, significant at the 3.5cr level, and 3.1er above ex- 
pected level of intrinsic anisotropy, at £ > 2010. Higher 
signal-to-noise ratio and multi-frequency measurements 
will be vital in confirming this signal and, if confirmed, 
determining its origin. The cosmological implications of 
the deep field results are discussed further in Papers V 
and VI. 

We thank Roger Blandford, Bryan Butler, Mark Gur- 
well, Marc Kamionkowski, Ken Kellermann, Rick Perley, 
Sterl Phinney, and Wal Sargent for useful discussions. 
We acknowledge the invaluable efforts of Russ Keeney, 
Steve Miller, Angel Otarola, Walter Schaal, and John 
Yamasaki at various stages of the project. We gratefully 
acknowledge the generous support of Maxine and Ronald 
Linde, Cecil and Sally Drinkward, Barbara and Stanley 
Rawn, Jr., and Fred Kavli, and the strong support of the 
provost and president of the California Institute of Tech- 
nology, the PMA division Chairman, the director of the 
Owens Valley Radio Observatory, and our colleagues in 
the PMA Division. This work was supported by the Na- 
tional Science Foundation under grants AST 94-13935, 
AST 98-02989, and AST 00-98734.The computing facil- 
ities at CITA were funded by the Canada Foundation 



14 



MASON ET AL. 



for Innovation. LB and JM acknowledge support by 
FONDECYT Grant 1010431. SC acknowledges support 
by CONICYT postdoctoral grant 3010037. We thank 



CONICYT for granting permission to operate within the 
Chanjnantor Scientific Preserve in Chile. 



APPENDIX 



NOISE ESTIMATION 

An accurate estimate of the noise in each visibility measurement is important for power-spectrum estimation. For 
the CBI we estimate the noise as follows. In each 16-min scan (8-min lead and 8-min trail) wc form the differences 
of corresponding lead and trail integrations. With an integration time of 8.4 s, and allowing for slew and calibration 
time, there are usually about n — 50 matched integrations in one scan. The mean differenced visibility for this scan 
Vj and rms noise Sj is estimated from the individual integrations in this scan Vy : 



i=l 

1 " 



■? = 

The mean and variance are estimated separately for real and imaginary parts of the visibility and the two variances, 
which should be equal if the instrument is working and correctly calibrated, are averaged together. In the complete 
dataset a visibility measurement at a particular (u, v) point is usually constructed from many (to) such scans taken 
under different conditions and with different baselines, so the noise may vary from scan to scan. If the noise on each 
scan were known a priori, the maximum likelihood estimator of the visibility could be formed by weighted average of 
the scans 



with weights Wj — l/s|, and this estimator would have a gaussian distribution with variance 

" 2 = eV (A2) 

However, this is not the case: when the weights are estimated from the data, they have their own sampling distribution 
and the distribution of estimator V is not gaussian. V (eq. [Al]) remains an unbiased estimator of the visibility, but 
equation (A2) gives a biased estimator of its variance. The bias depends on the values n and to and on the range of Sj 
from scan to scan; when data from different baselines are combined to form a single visibility estimate, the bias also 
depends on the relative correlator gains and antenna temperatures. We have chosen to use the estimator of equation 
(A2) for the variance, but to correct it for the bias. We have found that, for large m and equal Sj, equation (A2) 
underestimates the variance by a factor 

/ ~ 1 + _ ~ 1.04. 

n 

In simulations we find / = 1.06. We attribute the discrepancy to fluctuations in the noises and numbers of samples 
in the actual data, second-order corrections, and a known and understood 1% overestimate of the noises in our 
pipeline. The simulations, and the 6% noise correction which results from them, take all of these effects into account. 
We should therefore increase the variances computed using equation (A2) by 1.06. We actually used an an earlier, 
incorrect, estimate of 1.08 for the factor /, so we have slightly overestimated the noise (by 2%, comparable with the 2% 
uncertainty in the noise variance discussed in S I4.4.1J1 . This overcorrection will have caused a small underestimate of the 
CMB band powers. The cosmological parameter analysis of Paper V corrects for the effect of this noise misestimation, 
as will future analyses. 
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